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Motivated by a recent experiment on the non-equilibrium dynamics of interacting bosons in ladder-shaped 
optical lattices, we report exact calculations on the sweep dynamics of Bose-Hubbard systems in finite two-leg 
ladders. The sweep changes the energy bias between the legs linearly over a finite time. As in the experiment, 
we study the cases of [a] the bosons initially all in the lower-energy leg (ground state sweep) and [b] the 
bosons initially all in the higher-energy leg (inverse sweep). The approach to adiabaticity in the inverse sweep 
is intricate, as the transfer of bosons is non-monotonic as a function of both sweep time and intra-leg tunnel 
coupling. Our exact study provides explanations for these non-monotonicities based on features of the full 
spectrum, without appealing to concepts (e.g., gapless excitation spectrum) that are more appropriate for the 
thermodynamic limit. We also demonstrate and study Stiickelberg oscillations in the finite-size ladders. 



I. INTRODUCTION 

The unprecedented experimental tunability of ultra-cold 
atomic systems in optical potentials provides the opportunity 
to study non-equilibrium quantum dynamics in previously in- 
accessible regimes 0]|2]. Very recently, pioneering experi- 
ments [3-5 1 have started experimentally exploring issues re- 
lated to adiabaticity, a fundamental concept of quantum dy- 
namics [6[. Of particular interest is the quantitative character- 
ization of deviations from ideal adiabatic behavior when the 
sweep rate of a system parameter is small but not infinitesi- 
mal. Deviation from adiabaticity in many-body quantum sys- 
tems has also generated extensive theoretical interest JUG). 
This may be regarded as the slow-sweep counterpart of the 
much-studied quantum quench (infinite sweep rate) Q. 

In the experiment of Ref. Q, which is the motivation of 
the present work, adiabaticity is explored through slow ramps 
of the energy bias (relative potential energy) between two 
legs of a ladder. The ladder is formed by two coupled one- 
dimensional tubes or chains. Alternatively, one can think of 
the experimental setup as a ladder of dimers (double well po- 
tentials) coupled to each other, as in Figure [T^. The experi- 
mental sequence starts by loading all bosonic atoms in the left 
leg. The initially filled left leg can be lower in potential en- 
ergy than the right leg {ground state sweep) or higher {inverse 
sweep). Subsequently the potential energy is reversed linearly 
in time, Figure[lJ). In a truly adiabatic sweep, i.e., in the limit 
of large sweep time r, it is expected that all the bosons get 
transferred to the right leg. The fraction of particles in the 
right leg at the end of the sweep, nfl(r), the so-called transfer 
efficiency, characterizes adiabaticity. 

For the ground state (g.s.) sweep the transfer efficiency in- 
creases with sweep time r, as expected from the quantum adi- 
abatic theorem. In contrast, in the inverse sweep case, after an 
initial increase tir(t) decreases with r. Following Ref. (3), 
we call this non-monotonic behavior a breakdown of adia- 
baticity. Of course, for extremely slow sweeps (extremely 
large r), the adiabatic theorem assures us that Ufl(r) has to in- 
crease again to ti,r(t — > oo) = 1. For some r values, nn(r) 
also has a non-monotonic dependence on the intra-leg tun- 
nel coupling. The two non-monotonicities are shown through 
a selection of experimental data (provided by the authors of 
Ref. 0) in Figures [TJcd). 



In this work we analyze exactly the sweep dynamics in 
small ladder systems. We find that already the smallest sys- 
tem of two bosons in a two-rung ladder shows the two non- 
trivial non-monotonicities in transfer efficiency, with respect 
to sweep time and with respect to intra-leg tunneling constant. 
Analyzing the energy spectra of small ladder systems allows 
us to explain both effects, and gives us an alternate and useful 
perspective on the physics involved in sweep dynamics for 
the Bose-Hubbard ladder and the breakdown phenomenon. 
The sweep dynamics can be interpreted as a combination of 
Landau-Zener (LZ) processes [8] through a complex network 
of avoided level crossings. This exact study allows us to un- 
derstand the main phenomena without recourse to mean field 
treatments or low-energy effective descriptions. In addition, 
small-ladder systems show additional interesting features in 
the dynamics (Stiickelberg oscillations) which are averaged 
out in larger systems. 

In the presence of an overall trapping potential, even if the 
central ladders are long, there are off-center ladders which 
consist of a few rungs. Thus, although Ref. [3 1 focuses on long 
ladders, the results of the present study should directly de- 
scribe some part of the experimental measurements performed 
for Ref. (3). 

We will use the Bose-Hubbard Hamiltonian to describe the 
system (Figure[lJ): 

h=- U£ b UhR + J\\ E C w ) + hc - 
V i i ' tr / (i) 

i,a i 

where 6j )(7 , b\ are the bosonic operators for the site on rung 

i (i = 1 . . . L s ) and leg a = L, R, and n, iCT = b\ a bi M . The 
length of the ladder (number of rungs) is L s . The parame- 
ters Jj_ and J|| are the inter-leg and intra-leg tunnel couplings, 
respectively, U is the on-site interaction energy, and S(t) is 
the time-dependent energy bias between the two legs. We use 
J± = 1, measuring energy [time] in units of J± [H/J±\. For 
convenience H = 1, 

The initial state has all bosons on the left leg. In our com- 
putation we implement this by taking the initial state to be the 
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FIG. 1. (color online) (a) Each rung of our Bose-Hubbard ladder is a double-well potential (a Bose-Hubbard dimer). Parameters J±, Ju, U 
and S are indicated pictorially. (b) Ground state sweep and inverse sweep. The initial state has all bosons on the left site. (c,d) Representative 
experimental data from Ref. y), showing two non-monotonicities of the transfer efficiency tir(t). (c) Non-monotonic behavior of nn(r) 

function ("breakdown of adiabaticity"), for different axial z-lattice depths, (d) ur(t) at a constant value of r [ 7r ^ T = 0.53(3) j as a function 

of the axial ^-lattice depth (roughly: as a function of inverse J||), for different inter-tube hoppings. 



ground state of the Hamiltonian ([TJ with 5 = — 10 6 . Thus the 
initial state of the time evolution is not exactly an eigenstate of 
the Hamiltonian H (5 = ±<So) at t = 0, but for relatively large 
(5o (we typically use <5o between 20 and 40), the initial overlap 
with the relevant eigenstate of H(t = 0) is nearly unity. 

During a sweep the energy bias between left and right legs 
is changed linearly in time: 



26, 



6{0<t< t) 



zjp-t — Sq, ground state sweep 
-^Q-t + So, inverse sweep 



where r is a sweep time. The dynamics is characterized by 
the fraction of particles on the right leg 



n R (t) :=iV- 1 ^(7V fl (t)), 



(3) 



where N is the total number of bosons. In our calculations, 
the time evolution of the quantum state is obtained by numeri- 
cal integration of the Schrodinger equation and exact instanta- 
neous eigenvalues and eigenstates of the Hamiltonian ([TJ are 
calculated using full numerical diagonalization. 

We start in Section [IlJ with the sweep dynamics of a sin- 
gle Bose-Hubbard dimer, and show that the g.s. and inverse 
sweeps have different behaviors already for the dimer. Sec- 
tions |ljJJ and [IV] contain our main results on Bose-Hubbard 
ladders. Here we identify pertinent features of energy spectra 
of finite ladders. We present and explain non-monotonic be- 
haviors of the transfer efficiency tir(t), observed already in 
the smallest ladders. In Section [V] we explain oscillations of 
nnir) (Stiickelberg oscillations), which appear in finite-size 
ladders but get washed out in the thermodynamic limit. In 



Section VI we summarize results and mention open problems. 
In the appendices we describe the treatment of non-interacting 
Bose-Hubbard ladders and show some results for a larger sys- 
tem. 



II. BOSE-HUBBARD DIMER 

We begin with an L s = 1 ladder, i.e., a single Bose- 
Hubbard dimer. Spectral properties and dynamics of the dimer 
has been studied by various authors (e.g. ll9i fT7l '). Linear 
sweeps of the energy bias from — So to were studied in lfT6ll 
in the large U limit. Here, motivated by the expeiments (3), 
we focus on small to intermediate U and on the ±<5o — > T^o 
sweeps shown in Figure [TJ). We will show that the dimer al- 
ready shows marked difference between ground state sweeps 
and inverse sweeps, but there is no 'breakdown' phenomenon 
in the inverse sweep. 

Figure [2^ shows energy spectra for the dimer. There are 
N + 1 states. At finite U, the states at the top of the spectrum 
show a cluster of level crossings with tiny gaps, forming a 
characteristic and well-known swallow-tail shaped structure 

(e.g. EEMa). 

Figure |2j) shows the time evolution of the right site occu- 
pancy for ground state and inverse sweeps, for relatively large 
sweep time r = 40. In the non-interacting case ground state 
and inverse sweeps are identical, which is expected because 
upper and lower parts of the spectrum are equivalent. The 
adiabaticity is quite good for r = 40, as seen from nn(t) and 
the overlaps with eigenstates. 

For U > the ground state and inverse sweeps show dif- 
ferent behavior. For the g.s. sweep, the adiabaticity is even 
better at the same t than at U = 0, because the ground state 
is a bit more separated at U = 1. In contrast, for the same r, 
during the inverse sweep several states become populated by 
a sequence of Landau-Zener transitions, causing a decrease of 
boson transfer to the right site. To perform the adiabatic trans- 
fer, a sweep has to be carried out much more slower than for a 
ground state sweep. The reason for this difference is that the 
energy level structure at the top of the spectrum now contains 
many level crossings with small gaps. 

Oscillations in n R (t) (Figure [2] panels d) are due to occu- 
pancy of multiple eigenstates. There are already some oscilla- 
tions at the beginning of the sweep because the initial state is 
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FIG. 2. (color online) Spectrum and sweep dynamics of the Bose-Hubbard dimer (single-rung ladder) with N — 6 bosons, (a) Energy 
spectrum for U — and {7 = 1. (b,c,d) Sweep dynamics for (t,5o) = (40,20). For [7 = (first column), g.s. and inverse sweeps are identical. 
State numbering is shown for g.s sweeps and the inverse sweep simply has states 1,2,... replaced by 7,6,... (b) Population fraction in right site, 
nn(t). (c) overlap with the instantaneous eigenstates \ 4>i). (d) Frequency of oscillations of nn{t), contrasted with differences of instantaneous 
energy eigenvalues, (e) Transfer efficiency, showing no pronounced non-monotonicity or extended region of decrease. 



an eigenstate of H(S = -10 6 ) and not of H(S = ±S ). Dur- 
ing the sweep, eigenstates other than the lowest/highest get 
occupied further. Since the energy difference between instan- 
taneous eigenstates (and therefore oscillation frequencies) are 
changing with time during the sweep, we cannot obtain the 
frequencies by a Fourier transform of n#(i); in Figure [^pan- 
els d they are approximated from the time difference between 
successive minima. 

The transfer efficiency (Figure |2£) increases monoton- 
ically (neglecting oscillations) with t, for both types of 
sweep. Thus, in the dimer system inverse sweep dynamics 
shows no breakdown of adiabaticity, i.e., no pronounced non- 
monotonicity in tir(t). 



this simply means doubling the intra-leg tunneling constant 
J|| . Systems with periodic boundary conditions are transla- 
tion invariant and the dynamics considered in this paper pre- 
serves the total momentum along the leg direction. Since 
the initial state has zero total momentum, the dynamics is 
entirely confined to the zero-momentum sector. Thus we 
use only the subspace of the Hilbert space spanned by the 
linear combination of the position basis states of the form 

l/*/L^J2j=o 1 TJ \ n i,L,n 1}R , n Ls x, Ul^r), where T is 
the translation operator. 



A. Energy spectrum of Bose-Hubbard ladders 



in. BOSE-HUBBARD LADDERS 



In this section and in section IV we analyze the non- 
equilibrium sweep dynamics of ladder systems with a few 
rungs. In this section we describe features of the energy spec- 
trum, show what the sweeps mean in terms of the energy spec- 
trum, and demonstrate that the smallest possible ladder system 
(two bosons in a two-rung ladder, L$ — N = 2) already dis- 
plays the two intriguing non-monotonicities that are the main 
topic of this article. (Non-monotonic dependences of ur{t) 
as a function of r and as a function of Jy , closely analogous to 
the two non-monotonicities found in Ref. [3| and introduced 
in Figure [T]) In section IV we provide an explanation of the 
two non-monotonicities in terms of spectral features and tran- 
sitions between states of a diabatic basis. 

We will use periodic boundary conditions in the leg di- 
rection. For a ladder containing only two dimers (rungs), 



(Figure 
Figure 



Figure|3]shows energy spectra for N — 3 bosons in a L s — 
3 (three-rung) ladder. At large |<5 |, there are N + 1 groups 
of states ('bands'), corresponding to 0,1,2,. ..,iV bosons on the 
left leg. This is analogous to the dimer which has N + 1 
states. Within each of these bands, the states have splitting 
(dispersion) determined mainly by the tunneling within the 
legs, J|| . For U = J|| = 0, the bands are each degenerate 

to. 

3J3 shows that in the U — 0, Jn ^ case there are 
true (unavoided) level crossings (Figure|3jl inset). When both 
U and J11 nonzero, the crossings are all avoided (Figure 
inset). In the sweeps considered in this paper, going from z£ 
to =F^0j me system wavefunction crosses a complex network 
of avoided level crossings such as that shown in Figure [3]i 

In Figure [4^ we show the paths truly adiabatic sweeps 
would take, for both the g.s. sweep and the inverse sweep. 
In the two cases, the starting points are the lowest state of the 
lowest band (g.s. sweep) and the lowest state of the highest 
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FIG. 3. (color online) Energy spectra of Bose-Hubbard ladder for 
L s — N = 3 in zero momentum subspace. (a) U = Jii = 0: the 
eigenstates are degenerate, (b) Jm = 1, U = 0: each band splits into 
eigenstates with different kinetic energy. Level crossings are real (not 
avoided) crossings, as shown in inset, (c) U — 1, Jii = 0: the on-site 
interaction splits the eigenstates within each band, (d) Jm = U =1: 
all crossings in the spectrum are avoided, as shown in the inset. 
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verse sweeps are different (as for the dimer), and in addition 
there is also the breakdown phenomenon. In this subsection, 
we present the main features of the non-equilibrium dynamics 
in few-rung Bose-Hubbard ladders during linear sweeps of the 
energy bias, and also present the main features of the transfer 
efficiency. 

Figure |4j} shows the time dependence of the boson fraction 
on the right leg nji(t), for a r = 40 sweep. For this sweep 
time, the ground state sweep exhibits almost adiabatic behav- 
ior, while the inverse sweep deviates significantly. This dif- 
ference is similar to that described for the dimer (single-rung 
ladder) in section [II] Compared to the dimer, the fundamen- 
tally new feature in the ladder system is seen in the depen- 
dence of ^^(t) on the sweep time r, which is non-monotonic 
and shows a pronounced region of decrease with increasing t. 
This is displayed in Figure [5] 

Oscillations in nn(t) complicate the behavior of the trans- 
fer efficiency riij(r). So, after the sweep, we evolve the final 
state with the time independent final Hamiltonian H(t = r) 
over a large additional time, to get the time-averaged transfer 
efficiency nfl(r) := (nn(t > r)) t . 

Figure [5] shows tIr(t) as a function of the sweep time r 
and of the intra-leg tunnel coupling Jii for inverse sweeps in 
the smallest non-trivial ladder (L s = N = 2). The transfer 
efficiency exhibits non-monotonicities as function of both r 
and Jm. For faster sweeps n^r) increases with increasing 
r. In contrast, in the slow sweep regime (t > 20), tir(t) 
decreases with increasing r for Jn > 0.2. In addition, for 
r > 20, the transfer efficiency at fixed r also depends non- 
monotonically Jy. These two non-monotonicities are analogs 
of those observed in experiment (Figure [TJ, and are analyzed 
in detail in the next section. Compared to the larger systems 
considered in Ref. [4 |, we have additional oscillatory features, 
analyzed in section [V] 



FIG. 4. (color online) (a) Energy spectrum in zero momentum sub- 
space for L s = 3, N = 4, U = 1, Jy = 0.4, with 42 levels. 
The adiabatic paths for g.s. and inverse sweeps are highlighted, (b) 
Fraction of bosons on the right leg n_R (t) as function of time, during 
sweeps with (5o,t) = (20,40). For the g.s. sweep the nn(t) reaches 
nearly unity, following the adiabatic path closely. In contrast, after 
the inverse sweep tir (t) falls far short of unity, even with a sweep as 
slow as r = 40. 



band (inverse sweeps). 



B. Sweeps and adiabaticity in few-rung ladders 

As we found for the dimer, in Bose-Hubbard ladders inter- 
actions are necessary for seeing a difference in adiabaticity 
between ground state and inverse sweep. In the appendix [A] 
we provide analytical calculations for non-interacting (U = 0) 
ladders where the g.s. and inverse sweeps are identical and 
there is no breakdown phenomenon. 

In this subsection and the rest of the main text, we will con- 
sider finite interactions. For (7^0 ladders, the g.s. and in- 



IV. THE NON-MONOTONIC BEHAVIOR OF THE 
TRANSFER EFFICIENCY 

In this section we provide microscopic explanations of 
the two non-monotonicities of the transfer efficiency 7Tr(t) 
shown in Figure [5] The explanations are based on consid- 
erations of the structure of energy spectra of Bose-Hubbard 
ladders and the couplings between energy levels. The r non- 
monotonicity can be understood by neglecting the crossing 
structure in the small \S\ regime of the spectrum and only 
examining the coupling matrix elements between the large- 
\S\ states. This will be explained first in IV A The Jy non- 
monotonicity arises from features of the crossing structure of 



the spectrum and is described in IV B 



A. Non-monotonicity with sweep time 

In this subsection we describe the microscopic origin of the 
non-monotonic behavior of the n/?(r) function (breakdown 
of adiabaticity). We will detail an analysis for the simplest 
case of L s = N = 2 (two bosons in two-rung ladder, Figure 
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FIG. 5. (color online) Transfer efficiency nnij) as function of the 
sweep time r and the intra-leg hopping Ju for inverse sweep, for 
the "minimal" ladder L s = N = 2 and J7 = l,<5o = 30. Non- 
monotonicities with respect to both r and Jm are visible in the 3D 
plot, and are also shown separately in 2D panels. 



[6}, and extract the general explanation from this special case. 
In terms of the state numbering in Figure [6] adiabaticity in 
the inverse sweep involves starting at state 5 on the left and 
populating state 1 on the right. 'Breakdown' involves a non- 
monotonic dependence on r of the final state 1 population. 
We therefore examine and explain the overlaps of the final 
wavefunction with the final eigenstates, as a function of r. 

The details of the crossing structure of the spectrum at small 
| S | (ellipse in Figure|6^) vary widely with Hamiltonian param- 
eters (c.f. Figure [3]). On the other hand, the breakdown effect 
is robust through a large region of parameter space. This sug- 
gests that the phenomenon can be explained ignoring the exact 
crossing structure of the spectrum around 6 = 0. Our strategy 
is therefore to use the states at large 6 as a 'diabatic' basis and 
examine the coupling matrix among these states. 

Diabatic representation: The spectrum (Figure |6^) con- 
sists of three bands of two states each, and the initial state is 
very close to the lowest state of the highest band, i.e. state 5. 

The time evolution starts and ends at \6(t — {0, r})| ^S> 1, 
where the Hamiltonian ([TJ is dominated by the energy bias 
term <V2^i ~ If one regards the on-site inter- 

action (U) and the inter-leg hopping (J±) terms as couplings, 
then in analogy to the two-level Landau-Zener problem the 
eigenstates of the Hamiltonian with U = J± = can be re- 
garded as the diabatic basis. The spectrum of the coupling 
free Hamiltonian Hj ±= u =a is shown in Figure [6J3. The dia- 
batic basis can be calculated analytically using the procedure 
described in Appendix |A| for treating non-interacting Hamil- 
tonians. 

For large energy bias \5\ 1 the eigenstates of the full 
Hamiltonian coincide with the diabatic states. The coupling 
strengths between diabatic states sets the probability for tran- 




FIG. 6. (color online) (a) Spectrum of minimal ladder (L 3 = N = 
2), with Jm = 0.5, U — 1. The adiabatic path for inverse sweep (5 
— > 1) is highlighted with circles. The crossing structure for small |<5| 
is indicated by an ellipse. Since S(t = 0) > we use —5 on the hor- 
izontal axis so that the sweep goes from left to right, (b) The diabatic 
basis states, obtained by setting U = J± = 0. (c) Overlap of the 
final state with the eigenstates of the final Hamiltonian and transfer 
efficiency. Same parameters as in panel (a); 80 — 30. System goes 
into state 1 for fast sweeps, but population of state 1 is reduced with 
increasing r in favor of 'band-excited' states 4 and 6. (d) Larger in- 
teraction, for which there is no breakdown. The final population of 
state 1 increases monotonically (neglecting oscillations) with r. 



sition between these states. Loosely speaking, this means the 
overlap of the final wavefunction with a diabatic state after 
the sweep is larger when the coupling (either direct or higher- 
order) with the initial state is larger. In other words, a tran- 
sition with larger coupling should be noticeable already for 
faster sweeps (smaller r), while states weakly coupled would 
need a slower sweep to get populated. 

The coupling matrix between the diabatic states is obtained 
by representing the full Hamiltonian in the diabatic basis: 
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The diabatic basis states are ordered according to the number- 
ing of Figure [6] We note that states from neighboring bands 
with same kinetic energy along the ladder (i.e., states 1; 3; 5 
and states 2; 4; 6) are coupled by J±. States from the same 
band, if coupled, are coupled through the interaction U. We 
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will next show how these observations on coupling terms al- 
low an explanation of the breakdown phenomenon, and also 
allow us to predict regions of parameter space where the phe- 
nomenon does not occur. 

Breakdown of adiabaticity: In Figure[6}; we show the final 
overlaps as a function of r, for a parameter combination for 
which breakdown is exhibited. At fast ramp rates (small r), 
the states at the bottom of bands (states 1,3) get occupied, so 
that the occupancy of state 1 (and hence n^(r)) increases to a 
maximum. For slower sweeps, the occupancy of these states 
is lost to higher-momentum states in lower bands, (states 4,6). 
This decrease of state 1 occupancy (and hence of n#(r)) is 
the breakdown phenomenon. The same description holds for 
larger ladders (which have more bands and more states in each 
band) — (a) the band ground states are successively occupied 
at small times, so that the nij(r) and the lowest state of the 
highest band increase initially with r, and (b) for slower ramps 
excited states of lower bands gain weight at the expense of the 
highest band g.s. These processes are illustrated for a larger 
system in Figure [7] 

We can explain this in terms of two different coupling 
strengths (different time scales) for initial excitation of the 
highest band ground state (state 1) and for later transfer to 
excited states of lower bands. In the L s = N = 2 case, this 
relies on |{7 being smaller than v2Jj_, and for larger ladders 
we would require U to be smaller than c\ J± for some constant 
C\. The lowest states of successive bands can get populated 
with coupling strength Jj_, so that the time scale for initially 
populating the 'adiabatic' final state (state 1) is set by ■ 
The time scale for moving to excited states of lower bands 
is however set by the weaker coupling ~J 1 ^~ 1 U, and is thus 
slower, and can happen only for slower sweeps. This explains 
the loss of state 1 occupancy in favor of excited states of lower 
bands, at larger r. 

Population of band excited states is fundamental to the 
breakdown mechanism. This corresponds to the idea of lon- 
gitudinal low-lying excitations (gapless phonon excitations) 
along legs, highlighted in Ref. [3]. Note that our microscopic 
mechanism above does not rely on any idea of gapless exci- 
tations, which strictly speaking only occurs in infinitely long 
ladders. 

Parameter regimes without breakdown of adiabaticity: 

A strong confirmation of our picture is the absence of break- 
down at larger U, as shown in Figure [6jl. For larger U, the 
time scale for populating lower-band excited states is com- 
parable or smaller than the time scale for initial population 
of higher-band ground state. Not having a slower time scale, 
there is no longer any mechanism for Ur(t) decrease with r. 

Another regime with no breakdown is that of small Jn . A 
simple way to understand this is that, for Jn — > 0, the system 
corresponds to isolated single dimers, and we have seen in 
section|II]that the phenomenon is an absent in a single dimer. 



B. Non-monotonicity with intra-leg hopping 

To explain the non-monotonicity of jTr(t) with Jn at mod- 
erately large r, we now consider the crossing structure of the 




FIG. 7. (color online) Spectra for L s = N = 3, J\\ = U = 1, 
showing time evolution during inverse sweeps (<5o = 30) through 
overlaps, which are indicated by line widths. The width of lines is ~ 
ln(overlap), and if the overlap is smaller than 2% the line is dashed. 
Transfer efficiency has a maximum at r = 40, tir(t) = 0.87. As 
argued in text, we see that for small r (left) the main process is the 
successive excitation of band ground states. In the breakdown regime 
(right), the higher states of lower bands are excited. Arrows highlight 
this mechanism pictorially. 
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FIG. 8. (color online) Time evolution during inverse sweep in the 
minimal ladder (L s = N = 2), shown through overlaps indicated 
by line widths. The width of lines is oc overlap, and if the overlap is 
smaller than 5% the line is dashed. Here U = 1, r = 200, 5 = 20. 
Insets on right show blowups of the crossings in energy spectra, in 
the same scale, highlighting the differences between small and large 
Jii. 



energy spectrum at small \6\ and how this structure changes 
with Jj|. Figure [8] shows this for the minimal (L s = N = 2) 
ladder. 

For very small Jn (Figure [8h) the bands have little or no 
overlap, so the lowest state of the highest band has little or no 
proximity to any other energy level during the sweep. As a 
result, the transfer efficiency is large for moderate r. 

As Jn is increased, the width of each band increases, re- 
sulting in multiple crossings, in particular between the lowest 
state of the highest band and various excited states of lower 
bands (Figures [HJ? and [8};). Via these crossings, the wave- 
function can populate states of lower bands, leading to loss of 



transfer efficiency. 

When J|| is increased even further (Figure |HJl), the cross- 
ings happen at larger \8\ and have smaller gaps. (Physically, 
the gaps are smaller at larger \ S\ because the diabatic states are 
more sharply different and hence mix less.) Small-gap cross- 
ings behave more like true crossings. Thus at large Jy and 
moderate r the main path taken by the wavefunction ends up 
at the lowest level of the highest band, even though this path 
passes through many crossings, as shown in Figure[8}i. 

To summarize, increasing Jy has two effects: at first it in- 
creases the number of crossings, which reduces the transfer ef- 
ficiency, and then for even larger Jy it makes the crossing gaps 
smaller, which again increases the transfer efficiency. This ex- 
plains the overall minimum in the Tir(t) versus Jy behavior. 
Of course, the behavior is complicated by interference effects 
(Stiickelberg oscillations), prominent in the lower right panel 
of Figure B] In Figure [8] we have chosen four Jy values that 
highlight the overall trend of the 71r(t) versus Jy behavior. 

V. STUCKELBERG OSCILLATIONS 

In this section we analyze the oscillations observed in the 
transfer efficiency (Figures [5] [6] and |9^). While the overall 
trends discussed previously reflect the macroscopic physics 
which is the focus of Refs. (3]|4l, the oscillations are particular 
to our microscopic (small-system) ladders. 

When a system is swept through a parameter region con- 
taining multiple avoided level crossings, the wavefunction 
may split into two eigenstates at one crossing and, later in 
the sweep, levels with nonzero weight may meet again and 
interfere. Such interference can lead to oscillatory behav- 
ior of observables as a function of the sweep rate ll2D . The 
wavefunction component following the energy level Ei will 
accumulate the phase J Eidt. As a result, in quantities like 
Ur{t) there are interference terms carrying phase factors of 
type J [Ei — Ej] dt, for each pair of levels with nonzero 
weight. Since the sweeps change 6 linearly in time, dt — 
— 2j^"d<5. Thus the phase factors are 

JLf m S)-EM]*S = ^r (4) 

where A a 's are the areas enclosed between levels and cross- 
ings in the energy-^ plane, as shown in Figure [5J; inset. Thus 
the oscillations as a function of t carry frequency A a /25o- 

In Figure]^ the Fourier spectrum of Ur(t) is shown. There 
are three pronounced peaks. These frequencies can be com- 
pared to "enclosed areas" in the energy spectrum diagram 
(Figure |9j; inset). Figure |9j; shows excellent agreement be- 
tween the oscillation frequencies, obtained by Fourier analy- 
sis, and the exact numerically calculated areas A a multiplied 

VI. SUMMARY AND CONCLUSION 
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FIG. 9. (color online) Stiickelberg oscillations in minimal ladder 
{L s = N = 2), with S =20,U= 1, = 0.85. (a,b) Transfer ef- 
ficiency tLr(t) and its Fourier transform, (c) Frequencies compared 
to areas enclosed by crossing lines in the spectrum. Inset shows rel- 
evant region of spectrum and definitions of areas. 



arise in the approach to adiabaticity in interacting systems 
(3JIU. Motivated by this experiment, we have undertaken an 
exact study of Bose-Hubbard ladders of explicitly finite size. 
We have thoroughly examined the energy spectra, their cross- 
ings, and the transitions between them, and explained some 
non-trivial features of the transfer efficiency. 

The understanding emerging from this analysis is that, 
while interaction effects are vital for the "breakdown of adia- 
baticity" and related effects, the thermodynamic limit is not 
at all necessary for understanding the main effects qualita- 
tively. In particular, the two non-trivial non-monotonicities 
are both well-explained without appealing to concepts like 
gapless-ness of excitation spectra or nonlinearities in mean- 
field descriptions. 

The present analysis opens up various avenues of possible 
future research. Since a trapped system contains ladders of all 
possible sizes, one might ask whether from the data one can 
extract information relevant to the sizes we have considered 
here, in particular whether it is possible to study the interfer- 
ence features (Stiickelberg oscillations) which are prominent 
in small ladders. Second, non-equilibrium dynamics in the 
size regime between the present sizes and large sizes remains 
unexplored. This relates to the widely appreciated difficulty 
of interpolating between full quantum treatments of small sys- 
tems and mean-field treatments of larger systems. 

Appendix A: Non-interacting Bose-Hubbard ladder 

With periodic boundary conditions (translation invariance), 
the time dependence of the U = Bose-Hubbard ladder of 
arbitrary length can be reduced to a single 2x2 system, i.e., a 
single Landau-Zener type problem. 

Writing the Hamiltonian ([TJ with U = in momentum 
space, one can diagonalize using a unitary transformation for 
every momentum mode {q = 0, . . . , L s — 1): 



A recent experiment with trapped bosonic atoms in ladder- 
shaped optical lattices has highlighted the intricacies that can 



Otq = U q ■ a q x + V q ■ d q ,R 
Pq=V q - a q ^ L - U q ■ a q< R, 



(Al) 
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withu 2 +w 2 — 1. The diagonalizati on condition yields u/v 
2J±/5, and the resulting diagonal Hamiltonian is 



2tt 



H u=o = J2 -2J||Cos(^)(4a g + /3t/3 g ) 



(A2) 



The eigenstates can now be created by applying a|, /3| op- 
erators on the vacuum, e.g., the ground state at negative 5 is 
\ipg.s) = -^j (Po) N \ty- Dynamics is obtained using Heisen- 
berg equations of motion for the f3 q operator, which leads to 
the same 2x2 problem for each q mode: 



The Stiickelberg oscillations in tTr(t) are now far less pro- 
nounced, but otherwise the features of the breakdown phe- 
nomenon are very similar to those presented in detail for 



smaller systems in the main text. The overlaps of Figure 10 3 
show that band excited states, corresponding to longitudinal 
modes along the leg directions, are excited in the breakdown 



process, consistent with the picture developed in section IV A 



ih 





(A3) 



The fraction of bosons in the right leg is n^(t) — \u(t)\ 2 . 
The 2x2 dynamics is too simple for nji(t) to display any 
of the nontrivial behaviors we have studied in the main text 
with nonzero interactions. For large <5o and t, we can use 
the Landau-Zener formula to calculate the transfer efficiency: 
ti r (t) = 1 - exp[-27r J 2 /a], with a = ±2^. Of course, 
this increases monotonically with r. 




Appendix B: L B — 5, iV = 4 



In Figure 10 the breakdown of adiabaticity is shown for 
a larger system containing four bosons in five-dimer ladder. 



FIG. 10. (color online) (a) Breakdown of adiabaticity in larger 
system. L s — 5, N = 4, U = 1, do = 20. (b) Time evolution during 
an inverse sweep in the breakdown regime (r = 80. Jm = 0.5), 
presented (as in Figures[7]and[8]l through the line- widths representing 
overlap (~ In (overlap)). 
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